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O Abstract 

^3 We apply KAM theory to the equation of the forced relativistic 

pendulum to prove that all the solutions have bounded momentum. 
Subsequently, we detect the existence of quasiperiodic solutions in a 
generalized sense. This is achieved using a modified version of the 
Aubry-Mather theory for compositions of twist maps. 

O 1 Introduction 

00 

In this paper we are concerned with some aspects of the dynamics of the 
(T) differential equation 

| s(7T^p) +asm * = /(i) ' (1) 

where a > is a parameter and / : R — >■ R is a continuous and T-periodic 
real function satisfying 

f(t)dt = 0. (2) 

This equation, sometimes called the forced relativistic pendulum, has been 
considered by several authors. In [TH] Torres proved the existence of a T 
periodic solution after imposing some restrictions on the period and the size 
of /. Later, Brezis and Mawhin [3] proved the existence of a T-periodic 



solution for any /. The existence of a second T-periodic solution has been 
proved in (3j [2] . See also [TU E] for an alternative approach to the periodic 
problem. The equation (TO) can be seen as a relativistic counterpart of the 
classical Newtonian pendulum 

x + asinx = f(t). (3) 

This equation has been analysed from many points of view. In particular Levi 
[H] and You [20] proved that all the solutions of ^ have bounded velocity 
x(t) whenever pi) holds. The relativistic framework implies that \x(t)\ < 1 
and so the boundedness of the velocity is automatic. However we will prove 
that the results by Levi and You have a relativistic parallel when the velocity 
is replaced by the momentum 

x(t) 

pit) 



^l-i(t) 2 ' 

The main result of this paper says that if /(£) satisfies (J2J) then all solutions 
of ([TJ satisfy 

sup \p(t)\ < oo. 

Moreover we will prove that condition (J2J) is essential for this conclusion. In 

addition we will prove the existence of generalized quasi-periodic solutions 

with two frequencies 

2tt 
wi = — > w 2 e(-l,l) 

Note that we find solutions for each frequency u>2 and these solutions are 
quasi-periodic when the phase space of the pendulum is a cylinder. These 
solutions become subharmonic solutions when u\ and io^ are commensurable. 
To prove these results we consider the Hamiltonian formulation of (fTl) 
where the position q = x and the momentum p = if=^ are conjugate 
variables. After some changes of coordinates we will write the associated 
Poincare map in a form such that Moser Twist Theorem is applicable and 
so invariant curves exist. This property already implies the boundedness of 
the momentum. To apply Moser's theorem, estimates is some C k norm are 
needed. These estimates usually are tedious and cumbersome and one has 
to find the right way case by case. This is why trying to repeat the direct 
computation by Levi or the change of variable by You, one is lead to non 



trivial technical difficulties. Anyway, a more general technique, inspired by 
[18] and based on the differentiability with respect to the parameter of the 
solution of a differential equation, will simplify significantly the computa- 
tions. Moreover, it could also provide a simpler proof of the results of Levi 
and You for the Newtonian case. We also stress the fact that we can consider 
the case in which the period of the forcing and the period of the potential 
are not the same. Furthermore, the generality of this argument allows to 
consider a general nonlinearity g(x) in (flj. 

To prove the existence of periodic and generalized quasi-periodic solutions, 
one can use the theory of Aubry and Mather [T2j. In principle, to apply this 
theory we need to known that the Poincare map of equation (fTl) has twist. In 
the paper [XT] it was shown that it does not hold unless a restriction on the 
parameters is imposed, namely the condition a < |^ is necessary. Since we 
want to obtain results for arbitrary parameters we will apply a less standard 
version of Aubry-Mather theory. In [13] it is shown that the main conclusion 
of this theory still holds when the map is obtained as a finite composition 
of twist maps. The Poincare map II of equation (fij) can be seen as a finite 
composition II = f\ o • • • o f N where every /, is a " small-time- map" that is 
twist without any restriction. To consider finite composition of twist map is 
a great novelty in [13] and many results on twist maps admit an extension 
to this setting. In [T3] the twist of each map /, goes to infinity as the action 
goes to infinity. The relativistic effects prevents the velocity from being too 
large and this makes impossible to satisfy this assumption of large twist. For 
this reason we must modify Mather's theorem in order to adjust it to our 
situation. With this modified theorem we can produce periodic and quasi- 
periodic solutions whose oscillating properties are determined by the rotation 
number of the corresponding Mather set. 

2 Motions with bounded momentum 

Consider the equation 

it(v^w)- 9{x) = m (4) 

and assume that the functions / and g satisfy the following conditions 
(Al) g G C 7 (R), g(x + S) = g(x), J Q S g(x)dx = 



(A2) / e C(R), f(t + T) = fit), fif(t)dt = 0. 

where T and S are two positive numbers. Notice that when g(x) = — a sin a; 

and S = 2n we recover equation ([T|. 

Equation Q is in the Lagrangian framework. Actually it can be expressed 

in the form 

d fdL\ __ dL _ 

dt \dx ) dx 
where 



L(x, x, t) = -\/I -x 2 + G{x) + f{t)x. 

Here G represents a primitive of g. Notice that G is S'-periodic and of class 

C 8 . 

To our purposes, it will be convenient to pass to the Hamiltonian formulation, 



q = K 



p 



v^ (5) 

p^-H q = g(q) + f(t) 



with H(t,q,p) = a/I + p 2 — G(q) — f(t)q. We arrive to this system after 
having performed the classical Legendre transformation 

„ - g ( 6 ) 

From now on the symplectic coordinate p will be called the momentum. 
The Hamiltonian vector field (H p , —H g ) is bounded so all solutions of (TsT) 
are globally defined and the same holds for the solutions of Q undoing the 
change of variables. 

Notice that, due to the relativistic structure, the velocity of any solution is 
bounded and satisfies 

\x(t)\ < 1 for each iel. 

We will prove that also the momentum is bounded. This is equivalent to the 
more restrictive condition on the velocity, 

sup \x(t)\ < 1. 

teR 

Precisely 



Theorem 1. Assume that (Al) and (A2) hold. Then every solution (q(t),p(t)) 
of pD satisfies 

sup \p(t)\ < oo. 

Moreover, we will show that the null mean value of the function / is an 
essential condition in the above theorem. 

Proposition 1. Assume that (Al) holds and that f is a continuous and 
T -periodic function satisfying 

f-T 



f = fj f(t)dt?o. 



Then there exists R > such that if (q(t),p(t)) is a solution of |5p with 
\p(0)\ > R, the momentum satisfies 

lim \p(t)\ = oo. 

t—^oo 

Proof of these results will be presented in the following sections. More- 
over, we will perform the proof for the case 5=1, being conjugated to the 
general one through a change of scale. 

3 The approximated Poincare map 

The solution of (J5J) satisfying the initial condition 

q(0) = g , p(0) = Po 

will be denoted by (q(t; qo,Po),p(t; qo,Po))- The main tool of our work will be 
the Poincare map, the area preserving diffeomorphisms of the plane defined 
by 

IT : M 2 ^M 2 , U(q 0}Po ) = (q(T;q ,p ) } p(T;q 0} p )). 

The periodicity of g (remember that we suppose S = 1) implies that II 
satisfies 

n(g + l,Po) = n(g ,Po) + (l,0) 

and so IT induces a diffeomorphism of the cylinder TxK, where T = R/Z. 
On the other hand, the periodicity of / allows to describe the dynamics of 
system (p)]) in terms of the map IT. In particular, the condition 

sup|p(t;g ,Po)| < oo 



is equivalent to 



sup|p n | < oo 



nei 



where (q n ,Pn) = n n (g ,Po)- Similarly, 

lim \p(t;q ,Po)\ = oo 

t— >-oo 

is equivalent to 

lim \p n \ = oo. 

Notice that the boundedness of the vector field (H p , —H q ) plays a role in the 
proof of this equivalence. 

In view of this equivalences, to prove theorem [l] we shall look for non con- 
tractible invariant curves for the map II. Two disjoint invariant curves define 
an annulus that is invariant under the diffeomorphism II, so we can say that 
they act as barriers. Our aim will be to apply Moser's small twist theorem 
to the Poincare map II. With the promise of being more precise later on, we 
recall that Moser's theorem gives the existence of invariant curves for a class 
of maps of the cylinder whose lift has the form 

e^e+u + sfctf + Rifar)] m 

ri = r + dR 2 (6,r) {> 

supposing that the reminders R\ and R2 were small in some C k norm. Here 
5 plays the role of a small parameter. 

The coordinates (q,p) are not the best ones to have the Poincare map written 
in form rt7|), so perform the following symplectic change of variables 

q = Q 

p = P + G(q) + F(t) 

where F(t) is a primitive of /. Note that that F(t) is T-periodic and C 1 . 
We get the system 

■ q _ P+G(Q)+F(t) 

P = g(Q)(l- r I+Sm+£B ) w 

Now we can introduce the small parameter S > through the following 
change of scale 

Q = u, P=^ vE [1/2,7/2]. (9) 

6 



It is important to note that the strip Ex [1/2, 7/2] corresponds in the original 
variables to the time dependent region 

A & = {(q,p) el 2 : A + G(q) + F(t) < p < ? + G(g) + F(t)} 

and so from the boundedness of F and G 

p — >■ oo as 5 — > uniformly in v. (10) 

System (J8J) transforms into 

l+fa[G(u)+F(t)] 



U 



y /S^v 2 +(l+Sv[G(u)+F(t)]) 2 

-Sv 2 g(u)[l i+fa[GM+nt)] 



(ID 

x /52„2 + ( 1+(5 „[ G(u ) + _p( t )])2 J 

The change of variables (J9J) is not symplectic, but the Poincare map of sys- 



tems (11) is still conjugated to II. 



Note that if 6 = system (11 ) transforms into 



u = 1 
v = 

and taking any initial condition (uq,Vq) G R x (1/2,7/2) we have that the 
solution is well-defined for t e [0, T]. So, by continuous dependence, there 
exists A > such that the if 5 G [0, A] the solution is still well-defined for 
t G [0, T\. The coordinates (u, v) are the good ones to have the Poincare map 
written in form ([7]). To have a rough idea of why this is true, one can see 



trough a formal computation that system (11) has the following expansion 

for small 5 

u = 1 - \5 2 v 2 + 0(S 3 ) 

v = 0(5 3 ). 

Notice the fundamental fact that up to second order F and G do not play 

any role. Now one can obtain the Poincare map integrating and evaluating 

at t = T. 

We are going to make this argument rigorous and the key is the theory of 

differentiability with respect to the parameters. So, inspired by [18], let us 

recall some general facts. Consider a differential equation depending on a 

parameter 

! = *<«,.,*) (12) 

7 



where * : [0,T] x V x [0, A] ->■ M n is of class C°^+ 2 ^+ 2 7 1/ > 1 and V is 
an open connected subset of W 1 and A > 0. The general theory of differ- 
ential equations says that the solution z(t, Zq, 5) is of class C 0,u+2,u+2 in its 
three arguments. The following lemma will be crucial for our purpose, and 
generalizes the result [TBI Proposition 6.4]. 

Lemma 1. Let K be a compact set of T> such that for every zq G K and 
5 G [0, A] the solution is well defined in [0,T]. Then, for every (t,z,S) G 
[0, T] x K x [0, A] the following expansion holds 

dz 5 2 d 2 z 5 2 

Z(t, 2o, ^) = *(*> ZO, 0) + o" — (t, 2 , 0) + — — (t, 2 , 0) + -i?(t, 2 , 5) 

||i2(t,.,<J)||c, (x) -)>0 as 5^0 
uniformly in t G [0, T] . 

Proo/. For a function G C o '" +2 ' ,/+2 ([0,2 Y | x AT x [0, A]), the Taylor formula 
with remainder in integral form gives 

0(t, 2b, 5) = 0(£, 2 , 0) + ^|(t, 2 , 0)5 +-2Q^(t, z o, 0) + R 2 (t, zo, 5) 

where 

R 2 (t,z 0l 5) = 1 -J ^(t,z ,Z)(6-0 2 dt. 

Integrating by parts one gets 

R 2 (t,z ,6) = ±{2j 5 ^(t,z ,O(5-m-^(t,z ,0)5 2 } 

and through the change of variable £ = 5s we get 

i2 2 (Mo,<*) = 5 2 / (1 -s)[— (t,2o,5s) - — (t,z o ,0)]ds. 

from which it is easy to conclude using the regularity of the solution. □ 

Note that, by means of this lemma we have a semi-explicit formula for 



the solution of (12). This is very useful to compute its Poincare map. So, let 



us apply the previous lemma to system (11). First of all, calling Z = (u,v) 



system (11) can be written in the form 

Z = ^(t;Z,5). 

The initial condition will be denoted by Z(0) = zq = (uq, vq) and the cor- 
responding solution by z(t;zo,S) = (u(t;uo,Vo,S),v(t;uo,Vo,S)). We will 
suppose, by periodicity, that z E [0, 1] x [1,3]. From ^ we have that 

z(t;u ,v ,0) = (u + t,v ). (13) 



To compute the first derivative with respect to the parameter let us call 

dz i 



X(t; z , 5) = §f (£; z , 5). We need X(t; z , 0) that solves the Cauchy problem 



X = A(t)X + a(t) 
X(0) =0. 

where 

A(t) = -^(t- zit- z , 0), 0), a(t) = — (t; zit- z , 0), 0). 

A simple computation gives 

— (t;Z,0) = — (t;Z,0) = (14) 

so that 

X{t;u ,v ,0) = Q. (15) 

Now let us compute the second derivative. Let us call Y(t; z , 5) = |^f (t; z , 5) 
with components (£(£; Zo, 5),r)(t;z , 6)). We need Y(t;z o ,0) that solves the 
Cauchy problem 

(Y = A(t)Y + b{t) 

\Y(0) = 

where 

6 W = -qp& *(*; zoi °)> °) + 2 d^s^ z & *°> °)' °) x ( t; Zo ' °) 

+ ^(t; z{t- zo, 0), 0) [X{t; z , 0),X(t; z , 0)] 



and jy^(t;z(t;zo,0),0) is interpreted as a bilinear form from M 2 x M 2 into 



dz 2 ' 
I 2 . A simple computation gives 

d 2 q 



d5 2 



(t;z(t;z ,0),0) = (-v 2 ,0). 



From (14) and (15) we get the system 



i = -vl £(0) = 

7) = 0, 7/(0) = 



leading to 

r(t;no,^,0) = (-w 2 t,0). 

Next we apply lemma fl] using (13), (15) and (16). We have that 



(16) 



5 2 5 2 

Z(t; uo, v 0l 5) = (u + t, vq) + -ji-vfr, 0) + y#(*; «o, «o, 5), 

where the remainder R satisfies the estimate 

ll#(V,£)l|c5([o,i]x[i,3]) -> as 6-* 

uniformly in t e [0, T]. Finally, evaluating at t = T we get the following 
expression for the Poincare map 



Mi = m + T - y Tv o + t-RiKi, v , 5) 
vi = v + ^R 2 (u ,v ,S) 



(17) 



and 



\Ri{-, •,^)||c 5 (R/zx[i,3]) + H-R2G, •,^)||c8(K/zx[i ) 3]) -> as 5 -> 0. (18) 



4 Invariant curves vs. Lyapunov functions 

We saw that theorem [I] will be proved as soon as we could place any initial 
condition (qo,Po) between two invariant curves. In view of (10) it is sufficient 



to prove the existence of invariant curves for the map (17) as 5 — > 0. More 



precisely, we are going to prove the existence of a sequence of invariant curves 
T n approaching uniformly the top of the cylinder. Analogously one can prove 
the existence of a sequence of invariant curves approaching the bottom of the 



10 



cylinder. Finally we will prove proposition [T] to show that the null mean value 

of the forcing / is essential to have invariant curves. 

Concerning the boundedness, as anticipated, we have found the variables 



(u,v) in order to have the Poincare map written in form (17) and apply 
Moser's small twist theorem whose original version is in [TJ]. There are 
many versions of this theorem and we shall employ one coming from the 
works of Herman [El [7] and explicitly stated in [17]. To recall it, let T = K/Z 
and consider the infinite cylinder C = TxM and its strip A = T x [a, b] with 
b — a> | . The theorem deals with maps g : A — > C with lift 

9 1 = 6 + w + 8[a(r) + R 1 (6,r)] 

?"i = r + <5_R 2 (#,r) 

where a G C A [a, b], and Ri, R2 G C 4 (A). The number to G R is arbitrary and 
5 G (0, 1] is a parameter. Suppose that the function a satisfies 

Co 1 < a'{r) < c Vr E [a,b], ||a||c*[o,6] < c o 

for some constant cq > 1. Moreover, we suppose that g satisfies the intersec- 
tion property, in the sense that 

s(r)nrv0 

for every non-contractible Jordan curve T G A. 

Theorem 2 (|17|). Let g : A — >■ C be a mapping in the previous conditions. 
Then there exists e > 0, depending on only on c , such that if 

\\F\\ c i(A) + IICIIc 4 ^) < e 

the map g has an invariant curve. 

Now everything is ready for the proof of theorem [TJ Excepting for the 
intersection property, it is easy to see that the Poincare map expressed in 



the form (17)-(18) satisfies all the hypothesis of theorem 2 In this case, 



9 = u , r = v , u = T, a(v ) = — -^Vq and 8 is small enough. Concerning the 
intersection property, notice that from a result in [TTJ , the null mean value of 
/ implies that the Poincare map associated to system &fy is exact symplectic 
in the sense that the differential form 

Pidqi — pdq 
11 



is exact in the cylinder. Moreover, it is known that an exact symplectic map 
has the intersection property. Finally we can say that also map (17) has the 
intersection property because this property is preserved by conjugacy. So, 
an application of theorem [2] proves theorem [1} 

We have just proved that hypothesis (Al) and (A2) imply that the mo- 
mentum is bounded. To complete the study of the boundedness we need to 
prove proposition [TJ We will perform the proof supposing that 

1 <- T 



f=7 T J Q f(s)ds>0, 

the other case being similar. We just need to prove that there exists R 
sufficiently large such that if \po\ > R then the corresponding orbit of the 
Poincare map II is unbounded. In this case, a less subtle expansion of II, 
coming directly from system ([5]), will be sufficient. So, integrate ^ and get, 
for t G [0, T] 



q(t;q ,p ) =q + t + e(t,q ,p ) 

P(t] qo,Po) = Po + Jo g(q(s; q Q ,Po))ds + f* f(s)ds 



(19) 



where 



i , P(s]q ,p ) { 

e(t,qo,Po)= < . == -l)ds 

h [ Vl+P 2 {s;qo,Po) 

As pit] qo,Po) — > oo as po — > oo uniformly in q and t G [0, T], we have e — > 
as po — > oo, uniformly in q and t G [0,T]. 

Adding and subtracting J g(q + s)ds = G(q + t) — G(q ) in the second 
equation of ( 19 ) we get 

p{t; q , p ) =p + G(qo + t) - G{q ) + / f(s)ds + e(t, q , p ) 

Jo 

where 

z(t,qo,p )= / {g(qo + s + e(s,q ,Po)) - g(q + s)}ds. 
Jo 

The mean value theorem implies that e — > as po — > oo uniformly in q and 

t G [0, T]. Evaluating in t = T we get the following expansion of II: 

q x =qo + T + e(T,q ,Po) 

Pi=Po + G(q + T) - G(q ) + Tf + e(T, q , p ) 

12 



where e and e tends to zero uniformly in q as p tends to +00. 
Now, inspired by pQ, consider the function 

V(q,p) = q-G(p). 

and notice that 

V(U(q,p)) = V(q,p) + T(q,p) 

where 

r(q,p) = -G(q + T + e(T,q,p)) + G(q + T)+e(T,q,p)+Tf. 

Now, using the fact that G is bounded, one can find K such that if V(qo, po) > 
K then p is sufficiently large in order to have T(q ,p ) > -%-. For such a p 
we have 

V(U(q ,p )) >V(q ,p ) + J y >V*. 
So, by induction we can prove that 

Tf 

V(q n ,p n ) >V(q ,po) + n — , n>\. 

Finally we have that 

lim V{q n ,p n ) = +00 

n— >oo 

and remembering the definition of V and the boundedness of G we get that 

p n ->■ +00. 



5 Generalized quasi-periodic and periodic so- 
lutions 

We have just proved that all the solutions of ^ have bounded momentum 
and a natural question is to describe the kind of recurrent motions that can 
be expected. Periodic solutions of different types always exist ([H],[5]), and 
now we are going to look for quasi-periodic solutions. Precisely, we will prove 

Theorem 3. For every uj G (—T,T), there exists a family of solutions of 
J5p, X{(t) = (qz(t),pt(t)), with (6K such that 

X c+1 (t) = Xe(t) + (1, 0) and Xz{t + T)= X i+U {t). (20) 

13 



Moreover, the initial conditions 

£^g ? (0) and £ i-> p s (0) 
are 0/ bounded variation and 

lim *® = " 

t— >oo t T 

To understand why these solutions satisfy a kind of weak quasi-periodicity 
we define, inspired by [16J, 



It satisfies 

^(e t + t, e 2 ) = ^(0 1? e 2 ), ^(e 1? 2 + 1) = $ € (0i, 2 ) + (1, 0) 

and this says that the function $£ is doubly periodic once it takes values on 
the phase space Txl. The solution is recovered by the formula 

when $^ is continuous as a function of the three variables (£, 9i,02)- This 
function is quasi-periodic. Again we are assuming that it takes values on 
Txl. In the discontinuous case the solution will not be quasi-periodic in 
the classical sense but the bounded variation of the initial conditions implies 
that quasi-periodicity in the sense of Mather will appear. See [12J for more 
details. When the number u is rational, say u) = | with a and b relatively 
prime, then 

X i (t + bT)=X i (t) + (a,0) 

and the solution is periodic with period bT. Once more we are assuming that 
X% takes values on T x R. Classically these solutions are called subharmonic 
solutions of the second kind. 



To prove theorem |3j consider the change of variable 

Q = q 

P = p - F(t) 



14 



where F = f. System (ph transforms into 



P+F(t) 



y ^i+(P+nt)) 2 (21) 

P = g(Q) 



The Poincare map of system (21 ) has a particular form. Consider a partition 



of the interval [0, T] in N sub intervals of equal length 

L = — < , = (22) 

and consider the map Il LtT (Q , Pq) = (Q(t+L; t, Q , P ), P(r+L; r, Q , P )) = 



(Qi, Pi) where (Q(t; r, Q , P ), P(t; T,Qo, Pq)) is the solution of (21) with ini- 
tial condition (Qo, Pq) at time r. The Poincare map II of the system can be 
written as composition of such maps, precisely we have that 

n = n Ti0 = n Li(iV _i) L o • • • o n LiL o n Lj0 . 

So let us study such maps. It is worth recalling some definition inspired by 
[13]. Consider a C 2 diffeomorphism f(9,r) = (Q(8,r),R(9,r)) = (9i,ri) of 
the infinite cylinder TxM that is isotopic to the identity. Passing to the lift, 
the components satisfy 

0(0 + 1, r) = &(6, r) + 1, R{d + 1, r) = R(6, r). 

We stress the fact that in his work Mather required only a C 1 diffeomorphism, 
but for our purposes we will need more smoothness. The diffeomorphism is 
said 

• exact symplectic if the differential form RdQ — rdO is exact in T x M, 

• twist if dQ/dr > 0, while, if there exists /3 > such that dQ/dr > (5 
we will say that / is /3-twist, 

• to preserve the ends of the infinite cylinder, if R(8, r) — > ±oo as r — y 
±oo uniformly in 8, 

• to twist each end infinitely, if 0(9, r) —8 —y ±oo as r — y ±oo uniformly 
in 8. 

Now we can recall the 

15 



Definition 1. Let V 00 = {Jayo'Pp, where Vp is the class of C 2 diffeomor- 
phisms of the infinite cylinder that 

1. are isotonic to the identity 

2. are exact symplectic 

3. are /3 -twist 

4- preserve the ends of the infinite cylinder, 

5. twist each end infinitely. 

For our purposes we will need 

Definition 2. LetV p+,p, ~ be the class of C 2 diffeomorphisms of the infinite 
cylinder that satisfy properties 1., 2., 4. of the previous definition and 

3\ are twist 

5'. are such that 0(8, r) — 9 — y p± as r —y ±oo uniformly in 9, 

6. there exists M such that \R(9,r) — r\ < M for every (S,r)6TxK 
Now we can start the study of the map Hl, t - Notice that by the period- 



icity of (21 ) it can be seen as a map defined on the cylinder TxR. Moreover 



we have that 

Lemma 2. For every t e [0,T], the map Hl, t is exact symplectic in T x 

Proof. Inspired by [S] consider the function 

,■„-, P . i VL \ F 2 W + P(i;r,Qo,P )F(t) + l , r initio pii' h 



First of all, it follows from the periodicity of (21 ) that Q(t; r, Qq+1, Pq 



Q(t; r, Q , P ) + 1 and P(t; r, Q + 1, P ) = P(t; r, Qo, -Po)- Hence we have 

V T (Q + l,P ) = V T (Qo,Po) 
Now let us compute the differential dV T . 



16 



We have 



8Vr 

Wo 



T+L 



P 



dP 



[1 + (P + F) 2 } 3 / 2 dQ 
L ( p Qp 



g{Q) 



dQ 

dQo 



p dQ 



[1 + (P + F) 2 p/2 dQ dQ 



dt 



dt 



(23) 



first equation in (21) we get 

-T + L 



using the second equation in (21). Now, integrating by parts and using the 



P^dt 



dQ 



o 



[P 

[P 



9Q ]T+L 
dQ lr 

dQ 1T+L 



T+L P ^dt 



dQo 



0Q iT 
that, substituting in ([231 gives 



-+L 



P 



dP 



dQo 

Analogously we can get 
dV T 






L) 



[1 + (P + F) 2 ] 3 / 2 dQ, 



P{T) d-Qo {T) 



dQ 



P(r + L)-^(t + L) - P(t] 



dQ 



r). 



a 



dP - v ' J dPo y ' ' ' ' yJ dP 
Hence dV = P\dQi — PodQo and the lemma is proved. 

This is not the only property satisfied by the map. In fact we have 

Proposition 2. For every r G [0, T], we have H t ,l G V~ l,l 

Proof. First of all, from lemma |2] we have that the map Hl,t is exact sym- 
plectic and by a similar argument as in [11] condition (22) implies that for 
every r G [0, T], the map Hl, t is twist and isotopic to the identity. From 



equation (|21|) we have 

'Q{t;r,Q ,P Q ) 



Qo + £ 



P{s;t,Q ,Pq)+F{s) 
y/lHP(s;T,Qo,Po)+F(s)y 



--ds 



P(t; r, Q , P ) = P + £ g(Q(s; r, Q , P ))ds. 



Evaluating the second equation in t = t + L, the boundedness of g gives that 
H T; l preserves the end of the infinite cylinder. Moreover, evaluating the first 
equation in t = r + L and using the second we easily get 



lim (Qi 

Pq— 5>±oo 



Qo) 



±L 
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uniformly in Q . Finally, property 6. is a trivial consequence of the bound- 
edness of g. □ 



So, summing up we have that the Poincare map of system (21) can be 
written as a composition of maps in V~ L,L and this justifies the study that 
we are going to develop in the next section. 

6 Composition of twist maps and proof of 
Theorem 



Consider a finite family {fi}i=x t ...,N such that fi e V°° for every i. Let 
F — fx o • • • o f N . We have that F is a C 2 exact symplectic diffeomorphism 
of T x 1 that preserves the ends and such that twists the ends infinitely. 
However, it has not to be twist. 

In [13], Mather proved that one can associate to F a continuous function 
h(9,9i), called variational principle, that acts as a generating function for 
a twist diffeormorphism. The variational principle satisfies, among others, 
these relevant properties: 

(HI) h(9 + l,9 1 + l) = h{9,9 1 ), 

(H5) There exists a positive continuous function p on IR 2 such that 

h( 7 , 9,) + h(9, 71 ) - h(9, 9,) - hb, 7 i) > r r P 



if 9 < 7 and 9 1 < 71 , 

(H6a) there exists a > such that 

9 — > a9 2 /2 — h(9, 9i) is convex for every 9\ 
9\ — ¥ a9 2 /2 — h(9, #1) is convex for every 9. 

The function h in general is not differentiable but from (H6) one can prove 
that the one side partial derivatives dih(9±, 9±) and d2h(9, #i±) exist. Mather 
proved that there exist particular configurations (0j) that minimize an ac- 
tion. They are called minimal configurations and are such that the partial 
derivatives dxh(9i,9i + i) and c^/i^j-i, 0j) both exist and satisfy 

dxhtfi, 9 i+ i) + d 2 h(9^ 9i) = for every i. (24) 



This property allows to construct a complete orbit (0i,ri) of F defining 

fi = -dih(9~i,9 i+ i) = d 2 h{9i-x,9i). 

Once we have a minimal configuration (0j), we can define for (p, q) G Z x Z its 
translate T Ptq 9 by (T Piq 9)i = 9 i+q — p. In an analogous way we can define the 
translate of an orbit. The translate of a minimal configuration is minimal. 
Moreover, given two configurations 9 = (6*j) and T = (•y i ) we say that < T 
if 9i < 7j for every i. Two configurations G and F are comparable if either 
= T or > T or < T. Using these characterizations Mather proved 

Theorem 4 ([13]). Let F = fro- --o fa with fi G V 00 fori = 1, . . . ,JV. Then 
for every u G K t/iere exists an orbit (9i,fi) of F such that any two translates 
of (9i) are comparable and the sequence (9j) is increasing. Moreover, 

v di 

hm — = lo 

i— >oo % 

and ui is called rotation number. 

The connection between theorem [4] and the result in the first paper by 
Mather [12] is stated in the following 

Corollary 1. From the orbit (9i,fi) in the previous theorem, we can construct 
two functions (j), r\ : 1R — >■ M satisfying, for every t G M 

0(t + i) = 0(t) + i, v (t + i) = T]{t) 

F(<f>(t),T]{t)) = {<f>{t + u)Mt + v)) (25) 

where <ft is monotone (strictly if u ^ Q) and r\ is of bounded variation. 

Proof. Inspired by [15], let us consider, for every u, the set 

S = {iGR: t = ju - k for some (j, fc) G Z 2 }. (26) 

We have to distinguish whether u is rational or not. 

— If u is irrational, S is a dense additive subgroup of R and every pair (j, k) 

gives rise to a different number. We proceed by steps. 

STEP 1: definition of (j) on E. If t G £ we define 

0(t) = ^ - fc. (27) 
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We claim that the function <fi : £ — y R is strictly increasing: we have to prove 
that 

ju — k < j'oj — k' =>• Oj — k < Oji — k' 

that is, calling r = j' — j and s = k' — k, 

< ru) — s =>- Oj < 9j +r — s. 

The case r = is obvious, so suppose r ^0. Suppose by contradiction that 
for some j G Z 

Oj > j+r - s (28) 

we have, from the comparison property of the translated, that either 

0i > 0i +r — s for every i. 

or 

0i = 0i +r — s for every i. 

In the second case the orbit would be periodic and this is not compatible 



with an irrational rotation number. So, from (28) we can prove by induction 
that for every ngN 

Now suppose that r > 0. Taking the limit for n — y oo after having divided 

by nr we get 

s 
0>u- -. 
r 

that leads to a contradiction as we multiply by r. Notice that we can repeat 

the same argument and get the same contradiction for r < 0. 

Moreover, satisfies the periodicity property 

(f>(t + 1) = 4>(t) + 1 for each t E S. 

STEP 2: extension of <fi outside E. Given r£l-E, the limits 

6(t±) = lim (hit) 
v ' t-M-±,tes v ' 

exist and 0(r— ) < </>(r+). To extend to a monotone function on the 
whole real line it is sufficient to impose 0(r) G [0(r— ), «/>(r+)] and we choose 
0(r) = 4>(t—). In this way : R — )■ R is strictly increasing and satisfies 

<j>(t + 1) = (f>(t) + 1 for each t G R. 
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STEP 3: definition off] on E. Define, for t 6 E 

rj(t) = d 2 h((j)(t-uj),(j)(t)) (29) 

where fa is the variational principle associated to F. We claim that for every 
t,seE 

\rj(s) -r7(t)|<a(0(s)-0(t)) (30) 

where a comes from (H6a). Supposing t<swe have from the monotonicity 

4>(t — u)<(j)(s — u), <j){t) < (j)(s), (j)(t + u) < 4>(s + u). 

Inspired by [131 Proposition 2.6], we notice that if in (H5) we set 7 = (p(s—u), 
9 = (p(t — u), 6\ — 4>{t) — e, 71 = 4>{t) with e > 0, divide by e and let e — > 
we get 

d 2 h(<p(s - u), 0(t)-) < d 2 h{(j>{t - to), (f>(t)). 

remembering that the partial derivatives exist on the orbits. Moreover, from 
(H6a) and remembering that the one side partial derivatives of a convex 
function exist and are non decreasing, we have 

d 2 h(<p(s - u), <f>{s)) < d 2 h(<p(s - u), <P(t)-) + a(<P(s) - <P(t))- 

Combining these two inequalities we have 

V(s)<r ] (t) + a( ( p(s)- ( p(t)). 



Using using (24) we can see that also rj(t) = —dih((p(t), <p{t + a;)) so we can 
get analogously 

V(t)<ri(s)+<x{<f>(s)-<l>(t)) 

and conclude. 

STEP jj.: extension of rj outside E. If r ^ E we define 

T](t) = lim ri(i) (31) 

tfr.tes 



This is a correct definition. Indeed, from (30) we have that 

\r}{tn+k) ~ v{tn)\ < a\(p{t n+k ) - <p{t n )\, 
and, being (p{t n ) a Cauchy sequence, we have that r](t n ) converges and the 



limit (31 ) exists. In principle the limit could depend on the sequence. This is 
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not the case, indeed in case that ??(£*) — > h and r)(t%) — ¥ / 2 we can construct 
a new increasing sequence (r n ) having (£* ) and (t^) as sub-sequences. So also 
r)(T n ) has to converge to a limit that is the same as l\ and l 2 - So the definition 



(31) makes sense. With this definition we have that estimate (30) holds for 
every t, s £ R. Since is monotone and hence of bounded variation, we have 
that rj is of bounded variation. 
Now, from the periodicity property of h and (j> we get that rj(t + 1) = rj(t). 



STEP 5: property (25) holds. Let us assume first that t £ S. Then t = ju — k 
and 

0(t) = Bj - k, (j)(t + lo) = 9 j+1 - k. 

Moreover, 

n(t) = d 2 h{cj){t - w ), 0(t)) = a 2 /i(^_ 1; ^) = f,- 

and similarly i](t +u) — r~j+i- Since (9j — k,fj) is an orbit of F we conclude 
that 

F(<j ) (t),r ] (t)) = (<P(t + u),r ] (t + u)). 

Let us assume now that t £ R \ E. So we select a sequence (£ n ) converging 
to t with t n £ S and t n < t. Then we can pass to the limit in the identity 

F(<f)(t n ),r](t n )) = ((f)(t n + u),r](t n + uj)). 

The irrational case is done. 

— The case u = - rational is simpler. We can suppose that p and q are 
relative primes and that the corresponding sequence (0$) is periodic (in the 
sense that 0i +q = 9i + p). First of all notice that in this case, the subgroup 



S defined in (26) is discrete, precisely, 



£ = {-: deZ}. 
Q 

The representation t = ju — k is not unique, indeed t = j- — k = f- — k' 
whenever k' — k = Np and j' — j = Nq for some iV £ N. Anyway the 
periodicity of (9i) implies that 

j--k= f- -k' ^9 j -k = 9i> - k'. 

q q 



So we can define on E as in (27). As before one can prove that : E — > 



is increasing (non strictly). We extend it to a monotone function on the 
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whole Rasa piecewise constant function that is continuous from the left and 
taking only the values 9j — k. 

Finally, as before, one can prove that (pit + 1) = 4>(t) + 1. Moreover, 
the fact that <ft takes only values at points of a minimal orbit, we can define 
directly for tGl 

7 t (t) = d 2 h(<P(t-u),<P(t)). 



This function is of bounded variation and condition (25) is satisfied as well. 
To prove this we just have to repeat the same arguments as in the irrational 
case. Note that this time it is not necessary to pass to the limit. □ 

In our case, theorem [4] cannot be applied, as the hypothesis of the infinite 
twist at infinity is not satisfied. So we will present a modified version of the 
theorem. First we give the following notation: let T^ be a sequence of non- 
contractible Jordan curves that are invariant under a map /. This curves 
are called invariant curves. We say that T^ |" +oo uniformly if there exists a 
sequence r^ — > +oo as k — > +oo such that T^ C T x (r^, +oo). The reader 
can easily guess the meaning of T^ j. — oo uniformly. 

We can prove 

Theorem 5. Consider a finite family {/i}i=i,... ) jv where fa G V p+,p ~ . Let 
F = /i o • • • o f N . Suppose that F possesses a sequence (Tk) of invariant 
curves such that Tk t +°° uniformly as k — >■ +oo and Tk 4 — oo uniformly 
as k — >■ — oo Then, for every u> G (Np_,Np + ) there exist two functions 
0, f] : K — > R satisfying the same properties as in Corollary [IJ 

The proof of this theorem relies on the following lemmas 

Lemma 3. Consider f G V p+,p, ~ . Fix an interval [a,b]. Then, there exists 
f G V°° such that f = f onTx [a,b]. 

Proof. It is convenient to work with the generating function h(9, 9i). Remem- 
ber that it is a C 3 function defined on the set S = {p_ < 9\ — 9 < p + } C M 2 
such that h(9 + l,9i + 1) = h(9,9i) and satisfies the Legendre condition 
duh < 0. It generates / in the sense that the map / is defined implicitly by 
the equations 

'd l h(9,9 1 ) = -r 

d 2 h{6,e 1 )=r 1 . [ ' 

More details can be found in [8]. Notice that the strip T x [a, b] of the cylinder 
corresponds to the set £2 — { a (@) < 9\ — 9 < 0(9)} C S where a and (3 are 
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implicitly defined by 



-d x h(9,9 + a(9)) = a 
-8^(9,9 + 13(9)) = b. 



The functions a and are C 2 , 1-periodic and the Legendre condition implies 
that a(6) < (3(9). Moreover, we have that a(9) | p- as a — > — oo and 
(3(9) f p + as b — > +oo. Now take two larger strips Si = {a < 9\ — 9 < b} 
and S e = {a + e < 9\ — 9 < b — e} such that E2 C S e C Hi C S (cfr figure). 
Notice that, by compactness, there exists 5 > such that d^h < —5 on Ei. 




Now, fix e > small and extend d^h out of {p_ + e < 9^ — 9 < p + — e} as a 
C 1 bounded function (it is not important how you do it). So we can suppose 
that there exists a constant M\ > such that 



sup \d 12 h\ < Mi. 

(0 o ,0l)6K2 



(33) 
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Consider x a C°° cut-off function of IR 2 such that 

X = 1 on E e 

X = 0on{9 1 -9> b}. 

Moreover we can suppose that x = x(@i ~ @)i — X < 1 an d X > on 
{6 — e < 0i — < b}. Define the new function 

A = x di2h+( X -1)6. 

We notice that A G C\ A(0 X + 1,0 + 1) = A(0i, 0) and 

A = 9i2^ on E e 

A = -<5 on {9 1 - 9 > b} 

With a similar argument as in [10] we can consider the following Cauchy 
problem for the wave equation 

(d 12 u = A(9 1 9 1 ) 

I u(9,9 + h) = h(9,9 + ~a) 

{ (d 2 u - d x u){9, 9 + a) = (d 2 h - d x h){9, 9 + a). 

The solution h + is defined on the set {#i — 9 > a + e}, is such that h + G C 2 , 
h + (9 1 + 1,9 + 1) = /i + (#i,#), <9i 2 /* + = A and h + = h on E e . Now perform 
an analogous argument to modify <9i2/i also in the zone {9± — 9 < a} finding 
h~ . Finally glue h + and h~ through the common part E e to get a function 



h. Notice that di 2 h < —5 on IR 2 . The function h generates via (32) a 
diffeomorphism f(9,r) = {9\,r\) such that the relation 

001 1 



Or d 12 h 

holds. So the diffeomorphism / is /3-twist with /3 = l/max{— di 2 h} and sat- 
isfies property 5'. Moreover, as h = h on E £ , the diffeomorphism / coincides 
with /onTx [a,b]. D 

It is not hard to guess that we are going to use this lemma to modify 
the diffeomorphism F through its components /j. So, it is worth introducing 
some notation. Given / G V p ~ ,p+ and an interval [a, b] then the modified dif- 
feomorphism / with support [a, b] is the diffeomorphism coming from lemma 
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|3 Given F = f\ o ■ ■ ■ o f N with fa G V p ~' pjr , we will call F with support [a, 6] 
the diffeomorphism given by F = f\ o ■ ■ ■ o f N where every fi is supported in 
[a, b]. Moreover, notice that, if fi G V°° then trivially f\ = fi. Finally, F has 
coordinates (@(9,r),R(9,r)) while fi has coordinates (0^(6*, r), R^(9,r)) 
and the corresponding modifications have coordinates (0(0, r),R(9,r)) and 
(6«(0,r),£W(0,r)). 

Lemma 4. Consider f G V p ~' p+ . There exists K > snc/i t/iai /or ever?/ 
modified f with support [a, b] 

\R{0, r)-r\<K for every (9, r) e T x R 

uniformly in [a,b]. 

Proof. We have to prove that, given a modification with support [a, b], we 
have the estimate with the constant K independent on [a,b]. Consider the 
generating function h of /. We have to estimate the quantity 

\d 2 h{9,9 1 )+d 1 h{9,9 1 )\. 

Notice that, with the notation of the previous lemma, in [b — e, a + e] we have 
h = h so the estimate comes directly from property 6. in the definition of the 
class / e V p -' p+ . If 0i - 9 > b or 9 1 - 9 < a then R(6, r) = r and K = 0. So 
we only have to study the cases b — e<9\ — 9<b and a < 9i — 9 < a + e. Let 
us study the first, being the second similar. We need d'Alambert formula, 
valid for a function V G C 2 (1R 2 ): 

V(9,9i)=- ^ d v [ V d 12 V(Z,7i)d£ + V(9,9 + 5)+ 

Je+8 Je 

d 2 V(ri — 5, r))dr) 

e+s 

where 8 G R. Applying it to h and choosing 5 = b — e we get 

h(9,9 1 )=~[ 1 d v f A(Z,ri)dZ + h(9,9 + b-e)+ 

J9+b-e J9 



I d 2 h(r] -6 + e,rj)dr] 
Je+l-i 
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Let us compute the partial derivatives. The fundamental theorem of calculus 
gives 

d 1 h(9,e 1 )= / &(0 : r])dv + d 1 h(9,9 + b-e). 

Je+b-e 

Remembering the definition of A we have, integrating by parts 

-0i /-0i /-0i 

A(0, rj)drj = x (v ~ 0)d 12 h(9, r])drj + 5 / { X (v -9)-l}dr t = 

16+b-e Je+b-e J 0+b-t 

X (6i - 6)d 1 h{6, 0i) - d!h(9, 9 + b-e)- / x'(v - 9)d 1 h(9, rfidr] 

Je+b-e 

+ 5J 1 { X (V - 0) - l}d V 

Je+b-e 

where we used the fact that xip ~ e ) — 1- So 



dyh(p, e 1 ) = x (6i - e)d l h{e, ej + s { x (v -o)- iR 

Je+b-e 
-01 

x'(rj — 9)dih(9,r))drj. 

9+b-e 

Similarly, 

d 2 ~h(9, 9,) = x (9 1 - 9)d 2 h(9, 9 1 )-S f \x(0i ~0~ 1R 

Je 

+ \ , (0i-Od2h(^9 1 )d^. 

Now we can concentrate on the quantity 

\d 2 h(9,9 1 ) + dih(9,9 1 )\. 

To estimate it we first notice that 

\x(9 1 -9)d 2 h(9,9 l )+x(0 l -9)d 1 h(9,9 1 )\ = \ x (9 1 -9)\\d 2 h(9,9 1 )+d 1 h(9,9 l )\ < M 

using property 6 in the definition of the class V p+,p ~ . Moreover, with the 
change of variable 9\ — £ = 77 — 9 we get 



\8 T { x ( v -9)-l}dr ] -8 [' +C { X (9 1 - - 1R| = 0. 
Je+i-e Je 



-b+e 

[ X ( V - 9) - l}d V - 5 I 

9+b-e 
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So we just have to estimate the quantity 

'e+b-c 



[ ' he x{Qi - i)d 2 h{i, ejdt + / * x '(v - o)d 1 h{e, v )d v \ 

Je Je+b-e 



that, after the change of variable rj = C, + b — ein the first integral and having 
noticed that \x'\ is bounded, reduces to an estimate of 

\d 2 h{r)-b + e,ei) +8x11(9, r))\dr) 

+b-e 

< |0i-0-6 + e| max \d 2 h(r] - b + e, 9 X ) + 0^(9, r])\. 

6»+6-e<r/<6»i 

Remembering that we are working in the region b — e < 9\ — 9 < b, 

\9i-9-b + e\<e. (34) 

Now, by the Legendre condition, the function 

tf (77) = d 2 h(r) - b + e, 0i) + 5i/i(0, 77) 

is monotone, so max^ + j_ e< <01 |^/(^)| is either |^(0i)| or |\l/(0 + 6 — e)|. Sup- 
pose we are in the first case, being the other similar. We have 

|* (0i)| < \d 2 h(9 1 - b + e, 0i) - Sfe/i(0, 0i)| + \d x h{9, 00 + 9^(0, 0x)| 
< |9 12 /i(c,0 1 )||0 1 -0-6 + e|+M 



for some c £ [0,0i — b + e]. Now we can conclude using (34) and (33). □ 



Lemma 5. Let F(9, r) be a diffeomorphism o/Txt. Assume that F = f\ o 
• • • ° In with fi G V p+,p ~ for i = 1, . . . , N. Then, for every u> G (Np_,Np + ) 
there exists three non negative constant r* ; A and B such that 

6(0, r) — > u> + i] for r > r* 

6(0, r) — > a; + rj for r > r* 

6(0, r) — < w — i] for r < — r* 

6(0, r) — < w — 77 for r < — r* 

where F has support [— r* — A*,r* + B*] with A* > A and B*>B. 
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Proof. For simplicity of notation, let us prove it for N = 2. The proof goes 
by induction. If N = 1, then u G (p_,p + ) and then by property 5' in the 
definition of the class V p+,p ~ there exist r* > and rj > such that 

Q(9 , r) — 9 > u> + rj for r > r* 
Q(6,r) — 6 < u — i] for r < — r* 

Every modified F outside [— r*,r*\ is twist, so 

d(G(9,r)-9) 



Or 



>0 



and, remembering that F(9, ±r*) = -F(6*, ±r*) for every #, one can verify 

that also 

©(#, r) — 9 > u + r\ for r > r* 

Q(6,r) — 9 < tu — r] for r < — ?%,. 

Now suppose that F = fi o f 2 so that we fix u; G (2p_, 2p + ). From the case 
N = 1 there exist p* and ?? such that, for i = 1, 2, 



'e«(0,r)-0> ^ forr>p* 

e«(#,r)-#<^2 forr<-p, 
e^^r)-^^ for r<-p*. 



(35) 



Moreover, as f 2 preserves the end, there exists r* > p* such that R^ 2 \9, r) > 
p* for r > r*. So, for r > r* 

9(0, r) - = (1) (e (2) (#, r ), i? (2) (#, r)) - & {2) (9, r) + (2) (0, r) - > w + 7] 

Analogously we can suppose that 

0(9, r) — 9 < oj — rj for r < — r*. 

Now take the modified fi with support bigger than [— r* — K, r* + if] where 
if is the constant coming from lemma |4j Let us estimate the quantity 

0(#, r ) - 9 = (1) (e (2) (#, r), i? (2) (0, r)) - 9 (2) (0, r) + (2) (0, r) - 



for r > r*. It comes from (35) that 0^(0, r) — 9 > ^p-- It remains to prove 
that 

(l)(e(2) ( ^ r))j R(2) ( ^ r)) _@(2) ( ^ r) > <^7 
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If r* < r < r* + K then R^ (9, r) = R^ (9, r) > p* and we get the estimation 



through (35). If r > r* +K then, by the definition of K, we have R^ 2 \9, r) > 
r* > p* and we conclude as before. In an analogous way we have the others 
estimates. □ 

Lemma 6. Let F(9, r) be a diffeomorphism o/TxR that possesses an invari- 
ant curve T. Assume that F — f\ o ■ ■ ■ of N with f\ G V p+,p ~ for i = 1, . . . , N. 
Then, for every u> G (Np_,Np + ) there exist three non negative constants r*, 
A and B, such that the following holds. Let (9 n ,r n ) be an orbit of F or of 
a modified F with support [— r* — A*,r* + B*] with A* > A and B* > B. 
Suppose that for every 77 > we have 

9 9 

lim inf — < u + 77 and lim sup — > u — rj 

n^oo n n^oo n 

then there exists n G Z such that 

(On,r n ) G T x (-r*,r*). 

Proof. Let r* , A and B the constant coming from lemma |5j We can suppose 
that r* is large enough to have T C T x (— r*, r*). The invariant curve divides 
the cylinder in two components, the upper and the lower and both are F- 
invariant (resp. F- invariant). Notice that to prove it we must use the fact 
that F (resp. F) preserves the ends. It means that there cannot exist orbits 
that jump from the top to the bottom of the cylinder. So, if an orbit (9 n , r n ) 
of F or F is such that r n > r* for every n or r n < —r* for every n then 

lim inf — > u + 77 or lim sup — < u — i] 

n^oo n ' n^oo 77, 

respectively, in contradiction with the hypothesis. D 

Now we are ready for the 

Proof of theorem^ Fix u G (Np_,Np + ), consider the constants r*, A and 
B coming from lemma [6] By hypothesis, we can find two invariant curves 
T + and T_ contained, respectively in r > r* or r < r*. Let S be the compact 
region defined by such curves. Let F^' = /i ° ■ • • ° /j for j = 1, . . . , N. The 
sets F^\H) are compacts and so one can find a region E, defined by two 
invariant curves such that 

S U F (1) (S) U F (2) (S) U • • • U F (iV) (S) C intt. 
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Analogously, we can find and a region Si = T x [— r* — A*,r* + B*\ with 
A* > A and B* > B such that 

E U F«(E) U F^(E) U • • • U fW(E) C mffii. 

Now modify every /, outside the strip Ex applying lemma [3] and find the 
corresponding /j. So we get F = f\ o ■ ■ ■ o f n . The diffeomorphisms F 
satisfies the hypothesis of theorem UJ so we get an orbit (0 n ,r n ) of F with 

,fn) G E. Notice 



rotation number u. By lemma p^ there exists n such that 
that r + and T_ are also invariant curves for F and so by the invariance on 
E we have that (6 n ,r n ) G E for every n. But in E we have F = F so that 
(9m T n ) is also an orbit of F. Remembering corollary n] we get the thesis. D 

Finally, we are ready for 

Proof of theorem [3| From proposition[2]we can apply theorem[5]to the Poincare 
map II of system (21 ) and find for every u G (— T, T) two functions (f) and n 
such that 

<K£ + i) = 0(0 + 1, *K£ + i) =*?(£) (36) 

n(0(O,ry(e)) = (0(e + w), ^ + w)). (37) 



Let Xg(£) = (Q^(t),P^(t)) be the solution of (21) with initial condition 



(</>(£)) ^(O)- Notice that from (36) and uniqueness we have that 

X ?+1 (£)=X e (£) + (l,0) 



and from (37) and the definition of II, 

Xz(t + T)=Xz+„(t). 



so that (20) is verified. Finally, consider the limit 



lim 

t— >oo 



t 



We have that, for nT <t<{n + 1)T 

Qs(t) __ Q&) 



t 



t 



QjjnT) QsJnT) nT 
+ nT t 



where, being the vector field in (21) bounded, the quantity Q^(t) — Q^{nT) 
is bounded. So we can compute 



lim 

t— >oo 



Q&) 



lim 

n— >oo 



Qt(nT) 



lim 

n— >oo 



Qi+nu{fy 



t n->oo nT n->oo nT 

where [x] denotes the integer part of x and {x} 



lim r ^ +{ - }(0) + 



nT 
x — \x\. 



[nu\ 
nT 



UJ 

f 

a 
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